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Abstract. We present a microscopic model for nanoparticles, of the maghemite (7-Fe20a) type, and per- 
form classical Monte Carlo simulations of their magnetic properties. On account of Mossbauer spectroscopy 
and high-field magnetisation results, we consider a particle as composed of a core and a surface shell of con- 
stant thickness. The magnetic state in the particle is described by the anisotropic classical Dirac-Heisenberg 
model including exchange and dipolar interactions and bulk and surface anisotropy. We consider the case 
of ellipsoidal (or spherical) particles with free boundaries at the surface. Using a surface shell of constant 
thickness (~ 0.35 nm) we vary the particle size and study the effect of surface magnetic disorder on the 
thermal and spatial behaviors of the net magnetisation of the particle. We study the shift in the surface 
"critical region" for different surface-to-core ratios of the exchange coupling constants. It is also shown that 
the profile of the local magnetisation exhibits strong temperature dependence, and that surface anisotropy 
is reponsible for the non saturation of the magnetisation at low temperatures. 

PACS. 75.50.Tt Fine Particle Systems - 75.30.Pd Surface Magnetism - 75.10.Hk Classical Spin Models 



1 Introduction 

Surface effects in a nanoparticle are of great importance 
since they dominate the magnetic properties and become 
more important with decreasing size of the particle. In 
particular, the picture of a single-domain magnetic par- 
ticle where all spins are pointing into the same direction, 
leading to coherent relaxation processes, is no longer valid 
when one considers the effect of misaligned spins on the 

■ surface on the global magnetic properties of the particle. 
Indeed, even for strong exchange interactions the surface 

, is responsible for coordination defects, and in a particle of 
radius of the order of 4 nm, 50% (in the 7-ES2O3 nanopar- 
ticles) of atoms lie on the surface, and therefore the effect 
of the latter cannot be neglected. 

The magnetisation near the surface is generally lower 
than in the interior. This effect leads to a thermodynamic 
perturbation in exchange interactions near the surface, 
which can be sizable at high temperatures. In a finite 
system, the effect of temperature introduces, in addition 
to the usual renormalization of the spin-wave spectrum, 
a spatial dependence of the magnetisation (see [Q and 
references therein for the case of semi-infinite cubic crys- 
tals) . Indeed, the symmetry breaking at the surface results 
in a surface anisotropy. In most cases, this happens to 
be strong enough as to compensate for the work needed 
against the exchange energy that prefers full alignment, 
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and it is conceivable then that the magnetisation vector 
will point along the bulk easy axis in the core of the parti- 
cle, and will then gradually turn into a different direction 
when it approaches the surface. In addition to the ex- 
change interactions, which are the strongest interactions 
between atoms in a magnetic system, there are also the 
purely magnetic dipole interactions between the magnetic 
moments of the atoms and the interactions between the 
magnetic moments and the electric field of the crystal lat- 
tice (spin-orbit interactions) . The last two types of inter- 
actions are relativistic in origin and therefore correspond 
to energies that are much smaller than the exchange en- 
ergy on a short length scale, but as they are long range 
interactions they lead, in general, to non negligible contri- 
butions. Indeed, these interactions play an important role 
in that they introduce a preferred direction in the system 
as they correspond to a Hamiltonian (see below) that is 
not invariant under rotation operations. In other words, 
they lead to the appearance of an anisotropy energy, i.e. 
the dependence of the total energy of the system on the 
direction of magnetisation. Moreover, the most important 
feature of the relativistic interactions is that they change 
the magnetic moment of an atom from one site to another, 
and hence introduce a non uniform spatial distribution 
of the magnetic moment in the system. The treatment 
of the dipole-dipole interactions leads to two energy 
terms corresponding to the volume and surface charges. 
In a small magnetic system, such as a nanoparticle, only 
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the second contribution is important since it accounts for 
the shape anisotropy, and the former becomes negligible 
as one integrates over a small volume (see a detailed dis- 
cussion in H and references therein). 

Therefore, it is necessary to take into account all dif- 
ferent contributions to the energy, exchange and dipolar 
interactions, bulk and surface anisotropies, in any study 
of spatial distributions of the magnetisation in a small 
magnetic system. 

One of our ultimate goals is to understand the effect 
of surfaces on the thermodynamic and spatial behaviors 
of the magnetisation in small systems, since this is also of 
crucial importance to the study of the dynamics of such 
systems where the surface effects are considerable. In this 
case, the dynamics is rendered more complicated by the 
additional (metastable or stable) magnetic states corre- 
sponding to the surface configurations. 



Therefore, according to Mossbauer spectroscopy and 
high-field magnetisation measurements, we may think of a 
particle as containing a magnetically disordered surface of 
a certain width and a core more or less ordered depending 
on the size of the particle. The net magnetisation of the 
particle is given by the (weighed) sum of the contributions 
from the surface and the core. The increase of the surface 
contribution at low temperatures becomes more important 
with decreasing size of the particle, as shown in Fig.lc. 

In this work, we compute the thermal behavior of the 
different contributions to the magnetisation due to the 
surface and core for different values of exchange couplings 
and different surface-to-volume ratios of the number of 
spins. We also compute the spatial variation of the mag- 
netisation at different temperatures and the specific heat 
for different sizes. All calculations are performed in zero 
magnetic field. 



According to Mossbauer spectroscopy performed by 2 Model for magnetic particles 



Morrish et al. (see 
nal work by Coey 



for a review , and also the semi- 
the iron cations at the surface of 
the 7~Fe203 particles have a noncollinear magnetic struc- 
ture from 4.2 K to room temperature. The Mossbauer- 
spectroscopy analysis performed in || on 7-Fe 2 03 parti- 
cles show that the spectrum contains two components, one 
associated with the bulk and the other with iron atoms on 
the surface. The latter seems to disappear at temperatures 
in the range 30 — 75 K, depending on the mean diame- 
ter of the particle assembly. From these Mossbauer-effect 
analyses it was also inferred that the surface shell has a 
thickness of about 0.35 nm, independently of the particle 
volume. This fact will be used in our simulations to deter- 
mine the ratio of the surface and core number of atoms. In 
addition, measurements of the magnetisation at high fields 
performed on the 7-Fe2C>3 nanoparticles (see also || 
for cobalt particles) have shown that the magnetisation 
is strongly influenced by the surface effects, depending on 
the particle size. In Figs. [j] we show the field and ther- 
mal variation of the magnetisation for different particle 
sizes. In Fig. |l|a we see that there is a sudden increase of 
the magnetisation as a function of the applied field when 
the temperature reaches 70 K, and that the magnetisation 
does not saturate at the highest field value, i.e. 5.5 T. In 
Fig. [l]b there is an important increase of the magnetisation 
at low temperatures. In Fig.lc the thermal behaviour of 
the magnetisation at 5.5 T is shown for samples with dif- 
ferent mean diameter. There we see that the smaller is the 
mean diameter of the particle the more important is the 
increase of the magnetisation at very low temperatures. 
This can be explained by the fact that the surface compo- 
nent corresponds to a state with canted atomic moments 
at low temperatures [|| [||. Furthermore, open hystereses 
were observed Q in the nickel ferrite particles (NiFe2C>4) 
up to fields of 16 T, while the anisotropy field of the bulk 
sample is utmost 4 x 10~ 2 T. The authors consider that 
the core spins are colinear whereas those on the surface 
are misaligned. In the case of 7-Fe2C>3 particles, however, 
no such irreversibility in the magnetisation hysteresis loop 
was observed Ipl, |B|. 



2.1 Hamiltonian 

As was discussed in the introduction, it is well established 
§, D that any theory dealing with the spatial magnetisa- 
tion distribution must consider all three energy contribu- 
tions, exchange, anisotropy and magnetostatic, and that 
there is no obvious reason for neglecting any one of these. 
However, it should be emphasised that in very small par- 
ticles the exchange interactions are too strong to allow for 
subdivision into domains, so that the magnetostatic en- 
ergy term is considered here only in order to account for 
the shape anisotropy in the case of ellipsoidal particles. It 
should also be stressed, as discussed in the introduction, 
that because of the surface anisotropy and the symme- 
try breaking at the surface, there is no reason to consider 
the exchange interactions on the surface as equal to those 
in the core, and thereby the particle may acquire a non 
uniform distribution of the magnetisation. 

Therefore, to take account of all these energy con- 
tributions, we propose a model, based on the classical 
anisotropic Dirac-Heisenberg (exchange and dipolar) Hamil- 
tonian, describing the spinel structure of the (spherical or 
ellipsoidal) ferrimagnetic nanoparticles (e.g., 7-Fe2C>3). In 
this article the calculations are performed using classical 
Monte Carlo simulations of the magnetisation thermal be- 
havior of the particle with open boundaries at the surface. 

So, the exchange, anisotropy and Zeeman contribu- 
tions are given by the following anisotropic Dirac-Heisenberg 
Hamiltonian 

— , — , Nt 

n DH = -J2 E ^/3S«-Sf +n -^^(S 4 .e 4 ) 2 

i.n a,l3=A,B i=l 



(1) 



where J a p (positive or negative) are the exchange coupling 
constants between (the a, (3 — A, B) nearest neighbors 
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Fig. 1. a) Magnetisation as a function of the magnetic field of a diluted assembly of 7-Fe2C>3 nanoparticles with a mean 
diameter of 2.7 nm. b) Thermal variation of the magnetisation extracted from a) at a field of 55 kOe. c) Thermal variation of 
the magnetisation in a field of 55 kOe for three samples with different mean diameters (2.7,4.8,7.1 nm). 



spanned by the unit vector n; Sf is the (classical) spin 
vector of the a th atom at site i; H is the uniform field 
applied to all spins (of number N t ) in the particle, K > 
is the anisotropy constant and the single-site anisotropy 
axis (see definition below). A discussion of the core and 
surface anisotropy will be presented below. In the sequel 
the magnetic field will be set to zero. 

To the Dirac-Heisenberg Hamiltonian we add the pair- 
wise long-range dipolar interactions 

n, _ (9^b) 2 \ - ( s » ' s j) Rjj ~ 3 (Si ■ Rjj) ■ (Rjj ■ Sj) 
"■dip- 2 2^ ^5 

where g is the Lande factor, fj, B the Bohr magneton and 
Rjj the vector connecting any two spins on sites i and j 
of the particle, = ||Rij||- 

2.2 Method of simulation 

The particle we consider here is a spinel with two dif- 
ferent iron sites, a tetrahcdric Fe 3+ site (denoted by A) 
and an octahedric Fe 3+ site (denoted by B). The nearest 
neighbor exchange interactions are (in units of K) pf]| , || : 
J AB /k B =s -28.1, J BB /k B - -8.6, and J AA jk B ~ -21.0. 
These coupling constants are used in the Dirac-Heisenberg 
Hamiltonian H.dh in order to model the phase transition 
from the paramagnetic to ferrimagnetic order as the tem- 
perature is lowered down to zero through T c ~ 906 K. 
In the spinel structure an atom on site A has 12 nearest 
neighbors on the sublattice B and 4 on the sublattice A, 
and an atom on site B has 6 nearest neighbors on A and 
6 on the B sublattice; the number of B sites is twice that 
of sites A. The nominal value of the spin on sites A and 
B is 5/2, and this justifies the use of classical spins. We 
have also taken account of | of lacuna for each two B 



atoms randomly distributed in the particle. The nanopar- 
ticle we have studied contains N t spins (~ 10 3 — 10 5 ), and 
its radius is in the range 2-3.5 nm. Our model is based on 
the hypothesis that the particle is composed of a core of 
radius containing N c spins, and a surface shell surround- 
ing it that contains N s spins, so that N t = N c + N s . 
Thus varying the size of the particle while maintaining 
the thickness of the surface shell constant (~ 0.35 nm), is 
equivalent to varying the surface to total number of spins, 
N st = N s /N t , and this allows us to study the effect of 
surfaces of different contributions. All spins in the core 
and on the surface are identical, but interact via the, a 
priori, different couplings depending on their locus in the 
whole volume. We will consider both cases of identical in- 
teractions, and that of the general situation with different 
interactions on the surface and in the core. Although we 
treat only the crystallographically "ideal" surface, we do 
allow for perturbations in the exchange constants on the 
surface. This is meant to take into account, though in a 
somewhat phenomenological way, the possible defects on 
the surface, and the possible interactions between the par- 
ticles and the matrix in which they are embedded. In || it 
was assumed that the pairwise exchange interactions are 
of the same magnitude for the core and surface atoms, but 
there was postulated the existence of a fraction of missing 
bonds on the surface. On the other hand, we consider that 
the exchange interactions between the core and surface 
spins are the same as those inside the core. We also stress 
that we are only concerned with non interacting particles, 
so we ignore the effect of interparticle interactions on the 
exchange couplings at the surface of the particle. 

In our simulations we start with a regular box (X x 
Y x Z) of spins with the spinel structure having the prop- 
erties mentioned above, and then given the radius or total 
number of spins in the particle, we cut in a sphere or an el- 
lipsoid. We choose the center of the particle to lie at one of 
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where X a are random numbers satisfying —1 < X a < 
1. Then the energy change SE produced by this move is 
calculated from Eqs. (|l|) and li SE > 0, one compares 
exp(— SE/ksT) with a random number < ? < 1. If 
<r < exp(— SE/ksT) or SE < 0, the trial configuration is 
accepted as new configuration, otherwise it is abandoned 
and the old one is counted once more. It is shown that 
this prescription leads to 
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Fig. 2. Normalized distribution of the coordination number in 
a spherical and ellipsoidal particle for Nat = 40%. 



the lattice sites. In all cases free boundary conditions are 
used, and thus a spin is considered as a core or surface spin 
depending on whether it has or not its full coordination 
number. Fig. |^ shows the distribution of the coordination 
numbers in a spherical and ellipsoidal particle. Therefore, 
in our simulations each spin has a structure associated 
with it that contains the coordinates of the corresponding 
site, the nature of the site in the spinel structure (i.e. A 
or B), its locus (core or surface), the number of its near- 
est neighbours, and the total number of such sites. Thus, 
once the structure of each spin is defined, the spins outside 
of the sphere or the ellipsoid are discarded. This accom- 
plishes the simulation of the crystal structure and shape 
of the particle. Next, we proceed with the Monte Carlo 
calculations. 

The classical Monte Carlo method based on the Metropo- 
lis algorithm is now a standard method and detailed de- 
scriptions can be found in ; so here we summarize only 
the main procedure. We calculate the expectation value 
of a function /({S^}) of the spins Si 



(f) 



Tr [exp(-H/k B T)f] 
Tr aq>(-H/k B T) 



(3) 



by generating a Markov chain of spin configurations {S^} 
of the system and taking the average 



(4) 



From any initial (random in our case) configuration c we 
generate a trial configuration by choosing the spin coordi- 
nates of a randomly chosen lattice site i by (a — x, y, z) 



s° 



X a .A 



(5) 



(/> = lim /. 



(6) 



Since one uses a finite L, it is essential to choose the ar- 
bitrary step parameter A suitably so as to exclude an 
appropriate number of initial configurations from the av- 
erage (||), and to estimate the "statistical error" in a reli- 
able way. In our case, we always start with a configuration 
of randomly oriented spins and adjust A so that 75% of 
the moves were successful, and it turned out that about 
3000 Monte-Carlo steps per spin were sufficient to yield 
an accuracy of a few percent for the local magnetisation. 

It is well known, of course, that spontaneous symme- 
try breaking, which is usually accompanied by a change 
from a disordered state at high temperatures to a sponta- 
neously ordered state at temperatures below the critical 
one, can occur only in the thermodynamic limit. In a finite 
magnetic system the magnetisation at zero field 



I N 



(7) 



vanishes at all nonzero temperatures irrespective of the 
number N of spins in the system. In particular, even above 
the critical temperature (of the bulk system) a small mag- 
netisation is found by Monte Carlo calculations, due to 
fluctuations which in a finite system observed over a finite 
time have not completely averaged out |jff). However, at 
very low temperatures there exist, even in an inhomoge- 
neous system, clusters of spins aligned with respect to each 
other, and there should exist an intrinsic magnetisation. 
This is usually defined as follows 



M = 



\ 




(8) 



Of course, in a finite system the magnetisation M is nonzero 
at all temperatures. Even at temperatures well above the 
critical temperature one still obtains M cx 1/yN, and 
this leads to the appearance of a magnetisation tail at 
high temperatures. The magnetisation M tends to that of 
the bulk system (infinite size) in the limit N — > oo. In the 
present case, we calculate the magnetisation defined in (||) 
for the core and surface by summing over the correspond- 
ing spins. 

Anisotropy energies : In both cases of a spherical 
and ellispoidal particle, we consider a uniaxial anisotropy 
in the core and single-site anisotropy on the surface. The 
easy axis in the core is chosen along our z reference axis, 
and the sites on the boundary have uniaxial anisotropy, 



H. Kachkachi, A. Ezzir, M. Nogues,, E. Tronc: Surface effects in nanoparticles: application to maghemite 7-Fe2C>3 5 



with equal anisotropy constant K s , whose axes are cho- 
sen to point outward and normal to the surface (f|, fil|). 
More precisely, we define for each spin a unit gradient 
vector on which the spin magnetic moment has to be pro- 
jected. In the case of a spherical particle these anisotropy 
axes are along the radius joining the center of the parti- 
cle to the considered surface site. For an ellipsoidal par- 
ticle the easy direction in the core is taken along the 
major axis of the ellipsoid, which is also along the z di- 
rection. The bulk anisotropy constant was estimated by 
many authors (see e.g. p| and references therein) to be 
K\ ~ 4.7 x 10 4 erg/cm .In our simulations we normal- 
ized this constant to the number of sites and in units of K, 
k c = (K c /k B ) ^ 8.13 x 10~ 3 K, k B being the Boltzmann 
constant. On the other hand, the surface anisotropy has 
not been determined experimentally in magnetic oxides, 
but it has become clear, however, that the correspond- 
ing contribution is very large as compared with the bulk 
one. In our calculations we took k s = (Ks/ks) — 0.5 K; 
the corresponding surface anisotropy constant K s was es- 
timated to be ~ 0.06 erg/cm 2 (see pj, @). In @ k s 
was taken in the range 1 — 4 K, but for these higher val- 
ues of surface anisotropy the authors obtained high-field 
irreversibility whereas the experimental results (also pre- 
sented in H) show none for the 7-Fe2C>3 particles. 

Finally, in the case of a very small ellipsoidal parti- 
cle, as discussed at the beginning of sect. 2. A, the effect of 
the dipolar interactions TLdip boils down to a mere shape 
anisotropy, which is absent in a spherical particle. There- 
fore, the magnetostatic energy of an ellipsoid of semi-axes 



can be written as 



E — 



2V 



D x ■ Mi 



D z ■ Mi 



(9) 



where D a , a — x,y, z, are the demagnetizing factors, M a 
the components of the net magnetisation, and V the vol- 
ume of the particle, which is equal to N t in our calcu- 
lations. If all semi-axes are different it is not possible to 
express the D's in closed form. However, this is possible 
in the case of a prolate or oblate spheroid. In the case of a 
prolate spheroid, as is very common in permanent-magnet 
materials, i.e. Z > X = Y, the demagnetizing factors are 
given by JT3 



D z =4tt 
D r = D 



1 



2e 3 
1 



log 



1 



2e 



(10) 



y --(4n-D z ), 



where e = y 1 — (X/Z) , < e < 1, is the eccentricity of 
the ellipsoid. In these calculations we assumed that the 
easy axes of the magnetocrystalline and shape anisotropy 
are the same, though this is not the case in general. 

We have found that using the long-range dipolar in- 
teractions TLdip involving all possible pairs of atoms in the 
particle, or the macroscopic magnetostatic energy E mag 
yields within numerical errors, the same results, only that 
the former contribution is much more time consuming 
than the latter. 



2.3 Results and discussion 

2.3.1 Thermal variation of the magnetisation 

In Figs. |^a-c, we plot the computed thermal variation of 
the core and surface contributions to the magnetisation 
(per site) as a function of the reduced temperature r core = 
T/T c core , T c core being the highest core "critical tempera- 
ture" , for N t = 909, 2009, 3766 with N st = 53%, 46%, 41%, 
respectively. These values of N st have been determined by 
the fact that the thickness of the surface shell is constant 
(~ 0.35 nm), according to Mossbauer-effect analysis ||. 
They correspond to a diameter of circa 4, 4.6 and 6 nm, 
respectively. 

In Fig. ^d we plot the mean magnetisation defined as 
+ N c M core )/N t as a function of 
r core , for the same values of N st as before. The exchange 
coupling in the core (i.e. Jaa, Jab, Jbb), generically de- 
noted by J c , are taken as 10 times those on the surface, 
denoted by J s . Since there is no experimental estimation 
of the exchange couplings on the surface, the choice of 
such J s was guided by the fact that the critical temper- 
ature of the bulk material is circa 906 K, while the hy- 
pothetical "surface transition" occurs in the temperature 
range of 30 — 75 K (which is a factor of more than 10 
smaller) according to Mossbauer spectroscopy and high- 
field magnetisation measurements. Nevertheless, we have 
also considered the effect of other ratios J c / J 8 (see below) . 
We see that the surface "critical region" , corresponding to 
the would-be magnetic phase transition on the surface is 
in a range of temperatures lower than the critical tem- 
perature of the core. This is expected from the fact that 
since the molecular field acting on a surface spin is lower 
than the one acting on a core spin, as a consequence of 
the lower coordination number at the surface and hence 
the change in crystal field thereon M, [[l?], and also be- 
cause of the fact that J s are taken smaller than J c . This is 
in agreement with the results of Mossbauer spectroscopy 
and the magnetisation measurements at high fields, since 
the surface component has a "critical temperature" in the 
range 30 — 75 K while for the core component the crit- 
ical temperature is much higher (~ 906 K). Note that 
both transitions are smeared because of the finiteness of 
the system size (see discussion after Eq.(||)). We also note 
that the surface magnetisation M sur f ace decreases more 
rapidly than the core contribution M core as the tempera- 
ture increases, and has a positive curvature while that of 
M core is negative, as is most often the case. Moreover, it 
is seen that even the (normalized) core magnetisation per 
site does not reach its saturation value of 1 at very low 
temperatures, and this can be explained by the fact that 
the magnetic order in the core is disturbed by the relative 
disorder on the surface, or in other words, the magnetic 
disorder starts at the surface and gradually propagates 
into the core of the particle (see Fig. || below). The mag- 
netic disorder on the surface is, of course, enhanced by 
the single-site surface anisotropy which tends to orientate 
the spins normal to the surface. In Fig. ||d we see that 
the more important is the surface contribution the more 
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Fig. 3. Thermal variation of the surface and core magnetisation (per site) (Figs.3a-c) and mean magnetisation (Fig. 3d) as 
obtained from the Monte Carlo simulations of an ellipsoidal nanoparticle. The anisotropy constants are given in the text; the 
exchange interactions on the surface are taken to be 1/10 times those in the core. 



enhanced and rapid is the raising of the mean magnetisa- 
tion at low temperatures, and this behaviour bears some 
resemblance to Fig. [|c. 

In Fig. [| we plot the core magnetisation of an ellip- 
soidal nanoparticle with Nt = 909, 3766,6330 and the 
magnetisation of the isotropic system with the spinel struc- 
ture and periodic boundary conditions n as functions of 
the reduced temperature t pbc = T/Tj BC ', and N st = 
53%, 41%, 26%. Comparing the different curves, it is seen 
that both the critical temperature and the value of the 
magnetisation are dramatically reduced in the core of the 
particle. The reduction of the critical temperature is obvi- 
ously due to the finite-size and surface effects |ll[ . There 
is a size-dependent reduction of the critical temperature 
by up to 50% for the smallest particle. The same result 
has been found by Hendriksen et al.Jl^] for small clus- 
ters of various structures (bcc, fee, and disordered) using 
spin-wave theory. As to the magnetisation, the reduction 
shows that the core of the particle does not exhibit the 
same magnetic properties as the bulk material, and as 
discussed before, it is influenced by the misaligned spins 
on the surface. In Fig. [| we can also see that the higher 
N t the lower the magnetisation in the critical region and 
the higher the temperature at which the magnetisation 
approaches zero, and this is consistent with the fact that 
M oc 1/v^Vf at high temperatures, as discussed earlier. 
However, the increase of the critical temperature with 
Nt is not as clear-cut as it could be expected, and this 
can be understood by noting that the disordered surface 
(J s — J c /10, small coordination numbers, and single-site 
anisotropy) strongly influences the magnetic order in the 

1 This system is a perfectly ferrimagnetic material with peri- 
odic spinel structure and without vacancies, though such ma- 
terial does not exist in reality since all spinels present some 
degree of vacancy. This system will be referred to in the sequel 
as the PBC system. 
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Fig. 4. Thermal variation of the magnetisation of the PBC 
system with Nt = 40 3 , and the core magnetisation for Nt = 
909,3766,6330 with N st = 53%, 41%, 26%, respectively, as 
functions of t pbc (see text). The exchange interactions on 
the surface are taken equal to 1/10 times those in the core. 

core through the relatively strong exchange couplings be- 
tween surface and core spins, which are equal to those in 
the core. 



2.3.2 Effect of the ratio J c /J s 

In Fig. H we plot the thermal variation of the surface mag- 
netisation as a function of the reduced temperature T core 
for J s = J c /2, J c and 2 J c , and N t = 5269, N st = 40%. 
Here the thickness of the surface shell is greater than 0.35 
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Fig. 5. The computed thermal variation of the surface contri- 
bution to the net magnetisation of an ellipsoidal nanoparticle 
as a function of r core for J s — J c /2, J c , 2J C (see text). 



Fig. 6. Spatial variation of the local magnetisation of a spher- 
ical nanoparticle of 3140 spins, as a function of the radial dis- 
tance, for r core < 1, r core = 0.5, and r core ~ 1". 



nm, and is not taken from experiments, contrary to the 
values in Figs. |^. The reason for taking such N s t is merely 
to compare the results for surfaces with different thick- 
ness. At any rate, our aim here is to maintain N s t fixed 
to a given value and vary the ratio J c / J s . We see that the 
surface critical region is shifted to higher temperatures 
upon increasing J s , and only when J s = 2J C that both 
the core and surface magnetic "phase transitions" occur 
in the same temperature range, i.e. r core ~ 1. It is worth 
noting that the weaker the exchange interactions on the 
surface the lower the magnetisation of the latter. This re- 
sult remains the same upon lowering the surface width. In 
this case the number of spins having smaller coordination 
numbers, and hence weaker effective exchange energy, i.e. 
those spins on the outer shell of the particle, is very small 
as compared with the rest of spins in the particle. More- 
over, we may think that the spins on the outer shell follow 
the (strong) molecular field created by the other (inner) 
spins constituting a relatively ferrimagnetically ordered 
core. Setting J c = AJ S , we may determine the coefficient 
A at which the transition regions of the core and surface 
overlap. This may be done using the simple, though time 
consuming, cumulant method introduced a few years ago 
by Binder (see |ll|] for a review). 

2.3.3 Profile of the magnetisation 

We have determined the spatial variation of the local mag- 
netisation of a spherical particle of a fixed radius as we 
move from the center out onto the surface, at different 
temperatures. At nearly zero temperature (j core -C 1), 
the local magnetisation decreases with increasing radial 
distance in the particle. It starts from the saturation value 
of an iron atom (in fact, we computed the mean value of 
A and B atoms, to take into account the possibility of 



starting at either atom) at the center of the particle and 
decreases down to the value of a surface atom. Note that 
what is plotted in Fig. || is in fact the (normalized) projec- 
tion of the atomic magnetic moment along the easy axis 
(z-axis for the core, and normal for the surface spins), 
which is proportional to cos0, and thus Fig. || also shows 
the spatial evolution of the orientation of the magnetic 
moment inside the particle as the radial distance is var- 
ied. The decreasing of the local magnetisation confirms 
what was said before, that is even at very low tempera- 
ture the surface is in a magnetic order which is different 
from that in the core, and as discussed earlier, the mis- 
alignement of spins starts at the surface and gradually 
propagates into the core. This could be related with the 
gradual canting of spins confirmed by Mossbauer spec- 
troscopy [0, [|J, U even at 4.2 K. As the temperature 
increases (r core = 0.5,T core ~ 1~), the local magneti- 
sation exhibits a jump of temperature-dependent height, 
and continues to decrease. Since the local magnetisation 
depends on the direction of the radius vector, especially 
in an ellipsoidal particle, the curves in Fig. ^ present, in 
addition to the usual numerical inaccuracy (especially for 
T core ^ ]_-^ gome fluctuations due to the spatial non ho- 
mogeneities inherent to the lacunous spinel structure of 
the 7-Fe2C>3 particles for each two B atoms, randomly 
distributed). 

2.3.4 Specific heat 

In Fig. [?] we plot the computed specific heat of the PBC 
system (N t = 40 3 ) and that of an ellipsoidal nanoparticle 
with Nt — 909, 3766, 6330 (with the corresponding values 
of N s t given earlier). There is a transition marked by a 
sharp peak for the PBC system and a broad peak for the 
three particle sizes. This broadening is an obvious illustra- 
tion of the finite-size effects. Note also the shift by 50% to 
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Fig. 7. Thermal variation of the specific heat of the PBC 
system with N t = 40 3 , and that of an ellipsoidal nanoparticle 
with Nt = 909, 3766, 6330, and for the same energy parameters 
as in Figs. ^). 

lower temperatures of the critical region, as in the mag- 
netisation thermal behaviour in Fig. [|. We are intending 
to prepare appropriate samples for the specific heat mea- 
surements. To our knowledge, such measurements have 
never been performed on nanoparticle assemblies. 



3 Conclusion 

We have presented a microscopic model for magnetic nano- 
scale particles including exchange and dipolar interactions, 
and bulk and surface anisotropy, and investigated the ther- 
mal and spatial behaviours of the core and surface mag- 
netisation of a nanoparticle of different sizes, and hence 
of different surface contributions. We have found that the 
finiteness of the system size and free boundaries lead to a 
non uniform magnetisation profile decreasing toward the 
surface in the particle. Moreover, it turns out that surface 
anisotropy also leads to non saturation of the magnetisa- 
tion at low temperatures. 

In order to compare our results obtained for a single 
nanoparticle with our experiments on dilute assemblies of 
nanoparticles [Q, and in particular to compare with the 
curve of magnetisation vs. temperature given in Fig. |l|b, 
we have to include the Zeeman energy term in our Hamil- 
tonian and run our program for different particle sizes 
with randomly distributed easy axes. Unfortunately, this 
will require several days of CPU time, since a temperature 
sweep takes about 30 hours (without dipolar interactions) 
on a two-processor Alpha Work Station. It also remains to 
study the interplay between finite-size and surface effects 
in a nanoparticle of round shape (sphere or ellipsoid), as 
is the case here, since for a box-shaped particle one can 
treat separately finite-size and surface effects by consider- 
ing both cases of free boundaries where these two effects 



are mixed and that of periodic boundary conditions where 
only the former are present. In the present case, it can 
be seen from Fig. [l] that the magnetisation of the PBC 
system where only finite-size effects are present is larger 
than the bulk magnetisation, while the magnetisation of 
a (anisotropic) nanoparticle is lower. This implies, as was 
analytically shown in JTs| ] for a box-shaped particle with 
simple-cubic structure, that finite-size effects yield a posi- 
tive contribution to the magnetisation while surface effects 
render a larger and negative contribution, resulting in a 
net magnetisation that is lower than that of the bulk sys- 
tem. In [g(J it was also shown that the difference between 
the finite-size and surface contributions is enhanced by 
surface anisotropy, which leads to non saturation of the 
magnetisation at low temperatures (see Figs. ||||). Nev- 
ertheless, already from the present single-particle prelimi- 
nary study we can infer some conclusions about the effect 
of the magnetically disordered surface on the global mag- 
netic properties of nanoparticles. The surface contribution 
to the magnetisation presents a rather different behaviour 
from the core contribution. Furthermore, even at very low 
temperatures, the local magnetisation decreases with the 
distance from the center, showing that the magnetic state 
on the surface is definitely different from that in the core. 
We have also shown that there is a drastic reduction of 
both the critical temperature and the value of the mag- 
netisation in the core of the particle. The reduction of 
the critical temperature is obviously due to the finite-size 
and surface effects. As to the magnetisation, the reduc- 
tion shows that the magnetic properties of the particle 
core are different from those of bulk material, because of 
the strong influence of the misalignement of spins on the 
surface, driven by symmetry breaking of the cristal field 
and surface anisotropy, which propagates from the bound- 
ary into the core. Finally, we note that our results on the 
thermal and spatial behaviour of the magnetisation are 
in qualitative agreement with those obtained by Binder et 
al. |E1[ for spherical particles with simple-cubic crystalline 
structure. 
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